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ABSTRACT 

Tidal dwarf galaxies (TDGs) are recycled objects that form within the collisional debris of interacting/merging galaxies. 
They are expected to be devoid of non-baryonic dark matter, since they can form only from dissipative material 
ejected from the discs of the progenitor galaxies. We investigate the gas dynamics in a sample of six bona-fide TDGs 
around three interacting and post-interacting systems: NGC 4694, NGC 5291, and NGC 7252 (“Atoms for Peace”). 
For NGC 4694 and NGC 5291 we analyse existing Hi data from the Very Large Array (VLA), while for NGC 7252 
we present new Hi observations from the Jansky VLA together with long-slit and integral-field optical spectroscopy. 
For all six TDGs, the Hi emission can be described by rotating disc models. These Hi discs, however, have undergone 
less than a full rotation since the time of the interaction/merger event, raising the question of whether they are in 
dynamical equilibrium. Assuming that these discs are in equilibrium, the inferred dynamical masses are consistent with 
the observed baryonic masses, implying that TDGs are devoid of dark matter. This puts constraints on putative “dark 
discs” (either baryonic or non-baryonic) in the progenitor galaxies. Moreover, TDGs seem to systematically deviate 
from the baryonic Tully-Fisher relation. These results provide a challenging test for alternative theories like MOND. 

Key words, dark matter - galaxies: interactions - galaxies: dwarf - galaxies: irregular - galaxies: kinematics and 
dynamics - galaxies: individual: NGC 4694, NGC 5291, NGC 7252 


1. Introduction 

During the collision of two disc galaxies, tidal forces strip 
material out of the discs and form long tails of gas and stars 
(Toomre & Toomre 1972). Zwicky (1956) was the first to 
suggest that gas can concentrate in the tidal debris and col¬ 
lapse under self-gravity, leading to the formation of a new 
object: a tidal dwarf galaxy (TDG). Tidal debris around 
galaxy mergers often contains gas and stellar condensations 
that have masses, sizes, and star-formation rates (SFRs) 
comparable to those of dwarf galaxies (e.g. Weilbacher et al. 
2000; Boquien et al. 2010; Dabringhausen & Kroupa 2013). 
To be considered as hona-fide TDGs, these condensations 
must be self-gravitating entities that can survive for several 
Gyr against internal or external disruption (e.g. Due et al. 
2000; Bournaud & Due 2006; Recchi et al. 2007; Ploeckinger 
et al. 2014, 2015). For example, they should show a velocity 


* Based on observations made with ESO telescopes at 
Paranal Observatory under programme 65.0-0563, 67.B-0049, 
and 083.B-0647. 


gradient that is kinematically decoupled from the surround¬ 
ing tidal debris, pointing to a local gravitational potential 
well. Such kinematic decoupling has actually been observed 
in some TDG candidates, e.g. Bournaud et al. (2007, here¬ 
after B07) and Due et al. (2007, hereafter D07). Since TDGs 
are formed out of material pre-enriched in the progenitor 
galaxies, they are found to have higher metallicities than 
typical dwarfs (Due & Mirabel 1998; Weilbacher et al. 2003) 
and are often detected by GO observations (Braine et al. 
2001; Due et al. 2007; Boquien et al. 2011). The systematic 
deviation of young TDGs from the mass-metallicity rela¬ 
tion has also been used to identify older TDG candidates, 
which are no longer embedded in collisional debris (Hunter 
et al. 2000; Sweet et al. 2014; Due et al. 2014). 

A key aspect of TDGs is that they are expected to be 
free of non-baryonic dark matter (DM; Barnes & Hernquist 
1992; Elmegreen et al. 1993). This occurs because of two 
basic dynamical principles: (i) tidal forces have different ef¬ 
fects on the dynamically-cold disc and on the dynamically- 
hot halo, thus they effectively segregate baryons in the disc 
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Table 1: Galaxy Sample 


System 

RA 

(J2000) 

Dec 

(J2000) 

Vsys 

(kms“^) 

Dist 

(Mpc) 

s 

O 

Short Desorption 

Bona-fide TDGs 

NGC 4694 

12 48 15.1 

+10 59 01 

1160 

17 

1.5 

disturbed lenticular: post-merger 

VCC 2062 

NGC 5291 

13 47 24.5 

-30 24 25 

4378 

62 

14 

peculiar lenticular: head-on collision 

NGC 5291N, 5291S, 5291SW 

NGC 7252 

22 20 44.7 

-24 40 42 

4792 

66.5 

17 

advanced merger: spiral+spiral 

NGC 7252E, 7252NW 


(which end up forming tails, debris, and TDGs) from DM 
in the halo (which is too hot to form narrow tails); and 
(ii) once a TDG has formed, it has a shallow potential well 
with a typical escape velocity of a few tens of km s“^, thus 
it cannot accrete dynamically-hot DM particles with veloc¬ 
ity dispersions of ^200 km s“^. High-resolution numerical 
simulations, including gas and star formation, confirm that 
TDGs should be nearly devoid of non-baryonic DM (Due 
et ah 2004; Bournaud & Due 2006; Wetzstein et ah 2007). 

Few observational attempts have been made to mea¬ 
sure the DM content of TDGs. Braine et al. (2001) used 
CO linewidths to estimate the dynamical masses of eight 
TDGs and found no evidence of DM. These CO observa¬ 
tions do not resolve the gas distribution and kinematics, 
thus it remains unclear whether the CO linewidths can 
be interpreted as rotation in a local potential well. B07 
presented high-resolution HI observations of the collisional 
ring around NGC 5291, and identified three TDGs^ that are 
rotating and gravitationally bound. Assuming that these 
TDGs and the collisional ring have the same inclination 
angle (as inferred from numerical models), B07 measured 
dynamical masses a factor of ^3 higher than the baryonic 
masses, suggesting that TDGs may contain some uncon¬ 
ventional type of DM (either baryonic or non-baryonic). 
They hypothesized that the missing mass consists of very 
cold molecular gas, which is not traced by CO observations 
and is not actively forming stars, but that may constitute 
a large fraction of baryonic DM in galaxy discs (Pfenniger 
et al. 1994). On the other hand, Milgrom (2007) and Gentile 
et al. (2007a) showed that Modified Newtonian Dynamics 
(MOND) could reproduce the observed rotation velocities 
of these TDGs, providing an alternative explanation to the 
missing mass without any need of DM. As pointed out by 
Kroupa (2012), TDGs may help to discriminate between 
MOND and the standard model of cosmology, given that 
these two paradigms predict different internal dynamics for 
these recycled galaxies. 

In this paper we present a systematic study of a sam¬ 
ple of six bona-fide TDGs around three nearby systems: 
the late-stage merger NGC 7252 (“Atoms for Peace”), the 
collisional galaxy NGC 5291, and the disturbed lenticular 
NGC 4694 (see Fig. 1). Specifically, we build 3D disc mod¬ 
els to study the Hi kinematics of the TDGs and constrain 
their inclination angles. These techniques have been widely 
used and validated for typical dwarf galaxies (S waters 1999; 
Swaters et al. 2009; Gentile et al. 2007b, 2010; Belli et al. 
2012a,b, 2014b). For NGC 7252, we present new Hi data 
from the Karl G. Jansky Very Large Array (VLA) to¬ 
gether with long-slit and integral-field optical spectroscopy 
from the Very Large Telescope (VLT). For NGC 5291 and 
NGC 4694, we re-analyse existing Hi data from the VLA. 

^ Stricktly speaking, the recycled objects around NGC 5291 
are not tidal dwarf galaxies because they are not formed in tidal 
debris but in a collisional ring. For simplicity, we refer to them 
as TDGs since they are expected to share similar properties. 


In particular, we revisit the DM content estimated by B07 
for the TDGs around NGC 5291. 

2. Galaxy sample 

Table 1 lists the basic properties of the systems in our sam¬ 
ple and their bona-fide TDGs. All six TDGs satisfy the fol¬ 
lowing requirements: (i) they are “genuine” condensations 
of gas and young stars, which are not due to projection 
effects along the debris (see e.g. Bournaud et al. 2004); (ii) 
they have higher nebular metallicities than galaxies of sim¬ 
ilar mass, indicating that they are not pre-existing dwarf 
galaxies; (iii) they display a velocity gradient that appears 
kinematically decoupled from the tidal debris, pointing to 
a local potential well. For the TDGs around NGC 7252, the 
validity of these requirements is demonstrated in Sect. 4. 
For the TDGs around NGC 4694 and NGC 5291, we refer 
to D07 and B07, respectively. Here we provide a concise 
description of these systems. 

2.1. The disturbed lenticular NGC 4694 

NGC 4694 is a disturbed lenticular galaxy that inhabits the 
outskirts of the Virgo cluster. It shows asymmetric optical 
isophotes, dust lanes, and signs of recent star formation. 
A long Hi feature stretches to the West of the galaxy for 
^40 kpc and contains a TDG candidate: VCC 2062 (see 
Fig. 1). D07 concluded that VCC 2062 is a recycled galaxy 
based on the following observations: (i) a high metallicity 
(12+log(0/H):^8.6), (ii) a strong CO detection, and (iii) an 
Hi velocity gradient that is kinematically decoupled from 
the debris. Recently, deep optical images have been ob¬ 
tained as part of the Next Generation Virgo Cluster Survey 
(NGVS, Ferrarese et al. 2012), revealing that NGC 4694 has 
a strongly perturbed halo and a prominent tidal feature 
to the West (Lisenfeld et al. in prep.). This confirms that 
NGC 4694 is the result of an old merger (^0.5 Gyr ago). 
The Hi kinematics of VCC 2062 is modelled in Sect. 5.1. 

2.2. The collisional galaxy NGC 5291 

NGC 5291 is a peculiar early-type galaxy that lies in the 
outskirts of the cluster Abell 3574. It is surrounded by a 
giant Hi ring with a diameter of ^160 kpc and Hi mass of 
^5x10^^ Mq (Malphrus et al. 1997). A disturbed galaxy 
(the “Seashell”) lies at a projected distance of ~13 kpc (see 
Fig. 1), but is a fiy-by object with a higher systemic veloc¬ 
ity (by ^400 kms“^. Due & Mirabel 1998). According to 
the numerical model of B07, the HI ring was formed during 
a head-on collision ^360 Myr ago. Notably, B07 found that 
the “Seashell” galaxy is not massive enough to reproduce 
the Hi ring and identified the elliptical IC 4329 as a more 
likely interloper. The Hi ring hosts several TDG candi¬ 
dates with nearly solar metallicities (Due & Mirabel 1998). 
The three largest TDG candidates (NGC 5291N, 5291S, 
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Fig. 1: Optical images of the selected systems overlaid with FUV emission from GALEX (pink) and Hi emission from 
the VLA (blue). TDGs are indicated in yellow. North is up; East is left. The bar to the top-right corner corresponds to 
^10 kpc for the assumed distances (as given in Table 1). NGC 5291 {left): composite ESO-NTT {BVR) and SDSS image 
from B07. NGC 4694 (top-right): MegaCam/CEHT image obtained as part of NVGS (Eerrarese et al. 2012). NGC 7252 
(bottom-right): composite ESO-MPG (BR) image from Baade et al. (1999). 


5291SW) are associated with Hi velocity gradients that are 
kinematically decoupled from the underlying debris (B07). 
In Sect. 5.1 we model the Hi kinematics of these TDGs. 


2.3. The advanced merger NGC 7252 

NGC 7252 (“Atoms for Peace”) is a prototype late-stage 
merger from the Toomre sequence of interacting galaxies 
(Toomre 1977). Numerical models have been able to re¬ 
produce its morphology and kinematics by simulating the 
merger of two disc galaxies falling together ^600-700 Myr 
ago (Hibbard & Mihos 1995; Chien & Barnes 2010). The 
central remnant has a single nucleus and a de Vaucouleurs 
luminosity profile, suggesting the morphological transfor¬ 
mation into an elliptical galaxy (Schweizer 1982; Hibbard 
et al. 1994). Using VLA observations, Hibbard et al. (1994) 
identified two TDG candidates along the tidal tails, but 
the spatial and velocity resolutions of the Hi data were 
too low to study their internal kinematics. In this paper 
we show that (i) both TDG candidates are associated with 
a steep Hi velocity gradient that represents a kinematic 
discontinuity along the tail (Sect. 4.1), (ii) both galaxies 
have relatively high metallicities, confirming that they are 
not pre-existing dwarfs but are forming out of pre-enriched 
tidal material (Sect. 4.2), and (hi) the TDG to the North- 
West has a complex Ba kinematics, which is likely domi¬ 
nated by turbulent motions (Sect. 4.3). Eor both TDGs we 
present 3D disc models in Sect. 5.1. 


3. Observations and data analysis 

In this section we describe the collection, reduction, and 
analysis of new multiwavelength data of NGC 7252: Hi 
observations from the NRAO-VLA^ (Sect. 3.1) and op¬ 
tical spectroscopy using both EORSl multi-object spec¬ 
trograph and GIRAEEE integral-field unit (lEU) at ESO- 
VLT (Sect. 3.2). In Appendix A we compare the new Hi 
data with archival VLA data from Hibbard et al. (1994) 
and single-dish observations. Eor the existing VLA data of 
NGC 5291 and NGC 4694, we refer to B07 and D07, re¬ 
spectively. 

3.1. Hi observations 

Spectral-line Hi observations were carried out with multi¬ 
ple configurations of the VLA. We used the hybrid DnC, 
BnC, and BnA configurations to compensate for the fore¬ 
shortening due to the low declination of NGC 7252 (6 = 
—24? 7). Eor simplicity, we will refer to them as D, C, and 
B configurations. In order to achieve a velocity resolution of 
2.6 kms“^ across a velocity range of ^500 kms“^, we used 
the correlator in 2AD mode covering a 1.56 MHz-wide band 
with 128 spectral channels and tuning the two intermedi¬ 
ate frequencies (lEs) at different central velocities. One IE 


^ The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative 
agreement by Associated Universities, Inc. 
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Table 2: VLA observing parameters for NGC 7252 


Array configuration 

DnC 

CnB 

BnA 

Observations Date 

Jun 2008 

Feb 2008 

Sep-Oct 2007 

Bandwidth (MHz) 

1.56 (x2) 

1.56 (x2) 

1.56 (x2) 

Number of channels 

128 (x2) 

128 (x2) 

128 (x2) 

Time on source (hr) 

2.9 

5.0 

16.0 


was set at ^4825 kms“^, while the other one was set at 
^4570 kms“^. The two IFs were merged during the data 
reduction to form a single band centered at 4700 kms“^. 
The VLA observing parameters are given in Table 2. 

The t^i’-data were calibrated and combined using AIPS 
(Greisen 2003) and following standard procedures. Self- 
calibration was required for each configuration to re¬ 
move low-level artefacts due to a 1.16 Jy-bright source 
(PKS 2217-251) located to the South of the field at ( 0 ^ 2000 , 
^ 2000 ) = (22^ 20"^ 37.65^ -24°53'58.0"). This continuum 
source was subtracted using UVSUB, after combining the 
B+C+D data and merging the two IFs. The residual radio 
continuum was fitted with a first order polynomial consider¬ 
ing line-free channels, which were determined by inspecting 
channel maps based on the D-|-C data. The continuum was 
then removed from the B+C-hD data using UVLSF. Some 
faint residual emission is still observable in several chan¬ 
nels at the location of the bright source; this emission was 
systematically discarded during the masking of the channel 
maps (described below). The ?xr’-data were mapped using 
multi-scale cleaning (Cornwell 2008; Greisen et al. 2009) 
and different weighting schemes: (i) robust (RO) weighting 
(Briggs 1995) with the RO parameter set to 0.5, giving a 
beam size of 15.8" x 10.7"; and (ii) natural (NA) weighting, 
giving a beam size of 21.1" x 17.8". 

After the Fourier transform, the HI cubes were analysed 
using the Groningen Imaging Processing Systems (Gipsy, 
Vogelaar & Terlouw 2001). Total Hi maps were built by 
summing masked channel maps. After various trials, we 
used a mask obtained by smoothing the cubes to 25" in 
space and 12.5 kms“^ in velocity, and clipping at 2.5crs 
(where is the rms noise in the smoothed cube). Velocity 
fields were built from masked channel maps estimating an 
intensity-weighted mean velocity. Since the HI line profiles 
are complex and asymmetric, these velocity fields provide 
only an overall view of the Hi kinematics. Our kinematical 
analysis is based on full 3D models of the observations, as 
we describe in Sect. 5.1. 

3.2. Optical spectroscopy 

Optical spectroscopy was obtained as part of several ob¬ 
serving runs at ESO facilities. For both NGC 7252NW 
and NGC 7252E, low-resolution long-slit spectra (cover¬ 
ing the 4500-7500 A wavelength domain) were obtained 
in July 1994 with the red channel of the EMMI cam¬ 
era installed on the New Technology Telescope (NTT) at 
ESO/La Silla. Details on observing conditions and data re¬ 
duction are given in Due & Mirabel (1998). The spectrum 
of NGC 7252NW is shown in Eig. 2 (grey line). 

New observations were carried out in September 2000 
with the EORSl multi-object spectrograph installed on the 
Unit Telescope 1 (UTl) of the VLT. We simultaneously 
obtained the spectra of 19 objects, including NGC 7252E 



Eig. 2: Optical spectrum of NGC 7252NW obtained with 
VLT/FORSl (black, A < 5500 A) and NTT/EMMI (grey, 
A > 5500 A). The FORSl spectrum has been rescaled to 
match the H/3 line flux measured from the NTT spectrum. 

Table 3: VLT observing parameters for NGC 7252 

Date 22 June - 24 August 2009 

Instrument GIRAFFE /ARGUS 

Filter name HR15 

Grating Acentrai (nm) 665 

Spectral sampling (nm) 0.005 

Velocity resolution (kms“^) 2.3 

Field of view 4.2" x 6.6" 

Pixel size 0.3" x 0.3" 

Seeing ~0.65" 

Integration time (s) 55 x 1365 = 75075 


and NGC 7252NW. The position angle of the slit-lets was 
fixed to 72.9°. We used grism 600B (1.2 A per pixel) with 
a total exposure time of 4x1800 sec. The resulting spec¬ 
trum has higher spectral resolution and sensitivity than 
the EMMI one, but only covers the wavelength domain 
3700-5500 A. Data reduction was carried out using IRAE. 
The spectrum of NGC 7252NW is shown in Eig. 2 (black 
line). Spectrophotometric data and resulting metallicities 
are discussed in Sect. 4.2 and listed in Table 5. As is com¬ 
mon practice, the line intensities are normalized relative to 
H/3=100 and the errors are estimated using the IRAE splot 
procedure which relies on Poisson statistics. 

We also obtained lEU observations of NGC 7252NW 
to investigate its Ra kinematics. We used GIRAEEE on 
the UT2 of the VLT (Pasquini et al. 2002), obtaining 55 
exposures of ~23 minutes each. The observations are sum¬ 
marized in Table 3. The data were reduced using the ESO 
pipeline. Additionally, we corrected each reduced cube to 
heliocentric velocity, before combining all 55 exposures into 
a variance-weighted mean cube. This final cube was then 
inspected and analysed using Gipsy. To isolate the Ra emis¬ 
sion, we subtracted the continuum in the wavelength range 
6650—6680 A, avoiding [N ll] emission lines. The continuum 
was determined by fitting a first-order polynomial to line- 
free channels bordering the Ra line. A total Ra map was 
built by integrating the cube in velocity and clipping at 3a. 
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Fig. 3: Overview of the new Hi observations of NGC 7252. Top: total Hi map at 25" resolution overlaid on the ESO MPG 
i?-band image. Contours are at 0.3 (^3cr), 0.6, 1.2, and 2.4 Mq pc“^. The Hi beam is shown to the bottom-right. The 
bar to the bottom-left corresponds to 20 kpc. Middle: Hi velocity field at 25" resolution. The crosses mark the location 
of the putative TDGs. The dashed lines show the path followed to obtain the PV-diagrams. Bottom: PV-diagrams along 
the East tail {left) and Nort-West tail {right). The vertical lines correspond to the crosses in the velocity field. Contours 
range from 2cr to 8cr in steps of la = 0.7 mJy beam“^. See Sect. 4.1 for details. 
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Table 4: Properties of NGC 7252 and its tidal tails 


Component 

HI velocities 
(kms“^) 

Mhi 

(10® Mo) 

Mh2 

(10® Mo) 

Lb 

(10® Lo) 

Remnant 

4596-4618 

0.25 

3.5 

51.6 

ET 

4745-4835 

2.0 

- 

0.9 

NWT 

4660-4830 

3.2 

> 0.02 

2.9 

Global 

4596-4837 

5.45 

> 3.52 

55.4 


Notes. The Hi mass and velocities of the remnant refer to 
the so-called Western Loop (see Fig. A.l in Appendix A). 
Molecular masses are calculated assuming Aco = 2.0 x 

10^° (cm“^ / Kkms“^) and using the GO fluxes from Dupraz 
et al. (1990) and Braine et al. (2001). For the NWT the molec¬ 
ular mass is a lower limit because single-dish observations cover 
only a fraction of the tail. H-band luminosities are taken from 
Hibbard et al. (1994) and rescaled to a distance of 66.5 Mpc. 


4. New results for NGC 7252 

In the following we describe the results from the new multi¬ 
wavelength observations of NGC 7252. Specifically, we use 
(i) Hi data at 25" resolution to study the overall gas dis¬ 
tribution and kinematics, (ii) long-slit optical spectroscopy 
to determine the metallicity of several Hll regions along 
the tails, and (iii) IFU spectroscopy to investigate the Ra 
kinematics of the TDG to the North-West. 

4.1. Hi distribution and kinematics 

Fig. 3 (top panel) shows an i?-band image of NGC 7252 
overlaid with the total HI map at 25" resolution. In agree¬ 
ment with Hibbard et al. (1994), we find that the Hi gas 
is mostly located in the outer regions of the system. The 
Hi tails closely follow the stellar tails and are more ex¬ 
tended than the stellar tails in the available i?-band image. 
The East tail (ET) has a linear geometry with a twist in the 
outer parts, while the North-West tail (NWT) has a curved 
morphology with detached Hi clouds at its tip. In both 
tails, large Hi clumps match the location of the putative 
TDGs, which are identified as blue stellar condensations 
with Hll regions and UV emission (see Eig. 1). None of 
the TDG candidates are located at the tip of the Hi tails, 
thus they are not affected by projection effects but repre¬ 
sent genuine condensations of gas and stars (cf. Bournaud 
et al. 2004). Other Hi clumps are present in the tails but 
they are not associated with stellar over densities. 

The middle panel of Eig. 3 shows the Hi velocity 
field at 25" resolution. The ET and NWT display a reg¬ 
ular velocity gradient and are, respectively, red-shifted and 
blue-shifted with respect to the central remnant, having 
f^sys = 4740 ± 14 kms“^ (from CO observations; Dupraz 
et al. 1990). The velocity field in the NW part of the merger 
remnant is uncertain because two different kinematic com¬ 
ponents overlap at the same spatial position along the line 
of sight (see Eig. A.l in Appendix A): (i) the base of the 
NWT ranging from ^4700 to ^4850 kms and (ii) the 
so-called Western loop (Hibbard et al. 1994) ranging from 
~4500 to ^4650 kms“^. The masses of these different com¬ 
ponents are given in Table 4. We associate the Western loop 
to the remnant as it likely corresponds to tidal material that 
is falling back on the galaxy (Hibbard & Mihos 1995). 


The bottom panels of Eig. 3 show position-velocity (PV) 
diagrams obtained along the two tails (dashed lines on the 
velocity field). Remarkably, steep HI velocity gradients are 
observed at the location of the putative TDGs (correspond¬ 
ing to the vertical lines in the PV diagrams). In the NWT, 
the HI velocity gradient appears to be decoupled from the 
underlying tail, pointing to a distinct kinematical struc¬ 
ture. In the ET, the Hi velocity gradient follows the gen¬ 
eral trend of the tail, although it becomes slightly steeper 
at the location of the TDG. These anomalous kinematical 
structures are further investigated in Sect. 5 using data at 
higher spatial resolution and 3D disc models. 


4.2. Metallicities along the tidal tails 

Here we use the spectroscopic observations described in 
Sect. 3.2 to determine the nebular metallicity of the two 
putative TDGs and other Hll regions along the tidal tails. 
Direct estimates of the Oxygen abundance depend on the 
electron temperature of the ionised gas, which is probed 
by the [OIII] line at 4363 A. Unfortunately, this weak line 
could not be detected in any of our spectra, thus we esti¬ 
mate Oxygen abundances using semi-empirical methods. 

Eor the EMMI red spectra, the oxygen abundance 
can be estimated from the [Nil]/Ha flux ratio. We 
adopt the calibration of Denicolo et al. (2002) and find 
12+log(0/H) = 8.6 ± 0.1 for both TDGs. Eor the EORSl 
blue spectra, we use the R 23 method involving the [O ll], 
[O III] , and R/3 lines (after correcting all fluxes for internal 
extinction using the H/^/Hy decrement). We adopt the in¬ 
dependent R 23 calibrations of McGaugh (1991), as updated 
by Kuzio de Naray et al. (2004), and Pilyugin (2000). The 
R 23 method is degenerate and gives two values of the abun¬ 
dance for one observable. This degeneracy is broken by the 
[Nil]/[On] ratio (McGaugh 1994). This line-ratio is only 
available for NGC 7252NW and clearly indicates that the 
higher abundance is the correct one. The resulting abun¬ 
dance is remarkably consistent with that from the [N ll]/Ha 
method (see Table 5). The high metallicities confirm that 
both TDG candidates are not pre-existing dwarfs, but they 
are currently forming out of pre-enriched material ejected 
from the progenitor galaxies. 

Using additional spectra from the EORSl run, we esti¬ 
mate metallicities at two other locations: a giant Hll region 
near the tip of the ET and one Hll knot to the West of 
NGC 7252. We find Oxygen abundances ranging between 
8.65 and 8.8, as summarized in Table 5. The metallicities 
appear to be uniform along the tidal tails and close to the 
solar value, with no evidence of strong gradients over dis¬ 
tances of 80 kpc from the remnant. Eor comparison, the 
metallicity in the outskirts of isolated spiral galaxies may 
show gradients of several dex and reach 1/10 z© (Eerguson 
et al. 1998). This suggests that radial gas mixing occurs 
during galaxy mergers, in agreement with numerical simu¬ 
lations (Rupke et al. 2010a; Montuori et al. 2010) and recent 
observations (Kewley et al. 2010; Rupke et al. 2010b; but 
see also Werk et al. 2011). The latter works emphasise the 
gas dilution due to the transportation of metal-poor gas 
from the outskirts to the inner regions. Here, instead, we 
have evidence of an outwards gas transfer, where metal-rich 
gas has been brought to the outer regions. 
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Fig. 4: Ra distribution and kinematics in NGC 7252NW. Left: i?-band image (gray scale) overlaid with Hi (blue) and 
Ra (red) emission. Hi contours are the same as in Fig. 9 (top middle). The outer Ra contour corresponds to a lOcr 
detection. The cross shows the Hi dynamical centre determined in Sect. 5.1. Middle: enlarged view of the Ra emission. 
The circle to the bottom-left shows the Ra beam. Right: PV-diagram obtained from the Ra datacube along the two main 
peaks (dashed line in the middle panel). The vertical line corresponds to the circle in the middle panel. Horizontal lines 
indicate the Hi systemic velocity of NGC 7252NW (solid) and the velocity range covered by the Hi emission (dashed). 
The Ra velocity resolution is ^2.3 km s“^ (see Table 3 for details). 


Table 5: Spectrophotometric data and resulting metallicities in the tidal tails of NGC 7252. 



4730 


< 4710 

I 

4690 

'u 

o 

^ 4670 

JZ 

4650 
1 4630 

c 

4610 


Ha PV-diagram 


Region 

NGC 7252NW 

NGC 7252E 

Hii-ET 

Hii-W 

RA(J2000) 

22° 20"^ 33.6^ 

22° 20"" 56.P 

22° 20"" 59.3^ 

22° 20"" 41.2^ 

Dec(J2000) 

-24° 37' 24" 

-24° 4T 9" 

-24° 4T 11" 

-24° 40' 45" 

[OH]a3727 

331±66 


237±40 

256±49 

Rs 

28±5 


28±5 

24±4 

R^ 

47±5 


52±6 

47±5 

H/3 

100±2 

100±6 

100±5 

100±2 

[OHI]a4959 

37±2 

53±5 

41±5 

29±2 

[OHI]a5007 

109±4 

150±9 

125±5 

89±4 

[NII]a6548 

25±5 

23±6 



He 

285±30 

233±28 



[NII]a6584 

62±10 

58±11 



[SH]a6717 

51±10 

42±8 



[SH]a6717 

33±9 

33±9 



12-hlog(0/H): 

- [Nil]/He (Denicolo et al. 2002) 

8.6 

8.6 



- R23 (Pilyugin 2000) 

8.55 


8.65 

8.7 

- R23 (Kuzio de Naray et al. 2004) 

8.7 


8.8 

8.8 


4.3. Ha kinematics of NGC 7252NW 

For NGC 7252NW we also obtained IFU spectroscopy of 
the Ra line (Sect. 3.2). The GIRAFFE field of view covers 
only the innermost region of NGC 7252NW, corresponding 
to the peak of the Hi distribution (see Fig. 4, left panel). 
Since the VLT has a pointing accuracy of ^3", the astrom¬ 
etry was determined assuming that the Ra emission spa¬ 
tially coincides with the R-band emission, as suggested by 
their similar morphologies (see Fig. 4, middle panel). The 
Ra distribution is elongated towards the North-South di¬ 
rection and shows two peaks separated by ^3". 

The Ra emission is resolved in both space and velocity, 
but it does not show a coherent kinematic structure. Fig. 4 
(right panel) shows a PV-diagram obtained along the two 


Ra peaks (dashed line in Fig. 4, middle panel). The Ra line 
profiles are broad: the main Ra peak has a FWHM of ^80 
kms“^ (more than twice the Hi linewidth), while the aver¬ 
age FWHM over the entire field is ^45 kms“^. The central 
velocity of the main Ra peak (^4665 kms“^) is signficantly 
lower than the Hi systemic velocity of NGC 7252NW 
(4678±3 kms“^), but it is consistent with the Hi velocities 
on its approaching side. The Ra emission, indeed, does not 
coincide with the dynamical centre of NGC 7252NW but is 
shifted towards the Western, approaching side (see cross in 
Fig. 4, middle panel). In general, the ionized gas appears 
more asymmetric and disturbed than the Hi gas. Most 
likely, we are observing Hll regions dominated by turbu¬ 
lence and/or non-circular motions due to stellar feedback. 
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Table 6: Properties of the high-resolution Hi datacubes 


System 

Beam 

(asec X asec) 

PA 

n 

AF 

(kms“^) 

Rms Noise 
(mjy beam“^) 

Ref. 

NGC 4694 

14.7x14.3 

-70.65 

7.3 

0.6 

a 

NGC 5291 

8.0X8.0 

45.76 

10.4 

0.5 

b 

NGC 7252 

15.8x10.7 

62.96 

10.6 

0.6 

c 


Notes. To increase the signal-to-noise ratio, the robust cube 
of NGC 5291 (with a beam 6.6^'' x 4.36^') has been smoothed 
to 8" X 8", whereas the robust cube of NGG 7252 has been 
smoothed in velocity to ~10.6 kms“^. References: (a) Due et al. 
(2007); (b) Bournaud et al. (2007); (c) this work. 

5. Dynamics of Tidal Dwarf Galaxies 

In the following, we investigate the internal dynamics of the 
six TDGs in our sample using high-resolution HI cubes (see 
Table 6). These cubes are obtained using robust weighting 
during the Fourier transform and have higher resolutions 
than those obtained with natural weighting (Sect. 3.1), 
but at a slight cost in terms of signal-to-noise (S/N) ra¬ 
tio. Figures 5 to 10 provide an overview of each TDG: top 
panels show an optical image (left), a total Hi map (mid¬ 
dle), and a Hi velocity field (right); bottom panels show 
PV-diagrams along the major and minor axes from the ob¬ 
served cube (left), a model-cube (middle; see Sect. 5.1), and 
a residual-cube (right). The Hi maps were built by inte¬ 
grating the cube over a narrow velocity range (correspond¬ 
ing to the kinematically decoupled velocity gradients) and 
clipping at Icr. The Hi velocity fields were built by estimat¬ 
ing an intensity-weighted-mean (IWM) velocity considering 
the same narrow spectral range, and clipping at 3cr. 

Using the high-resolution data, all six TDGs are spa¬ 
tially resolved with at least 2-3 beams along the major axis. 
The diffuse HI emission associated with the underlying de¬ 
bris is resolved out by the VLA at these high spatial resolu¬ 
tions, but some Hi clumps are still visible near the TDGs. 
We warn that the HI maps and velocity fields are uncertain 
due to the low S/N ratio, low number of beams across the 
galaxy, and possible contamination from the underlying HI 
debris. Despite these limitations, the Hi morphology and 
kinematics resemble a rotating disc. These putative discs 
are studied in detail in the next section using 3D kinemat- 
ical models. 

5.1. Kinematical models 

To test the hypothesis that TDGs are associated with an HI 
disc, we built 3D kinematical models following similar pro¬ 
cedures as Lelli et al. (2012a,b, 2014b). As a first step, we 
inspected channel maps and PV-diagrams to obtain initial 
estimates of the dynamical centre (xq, yo)^ position angle 
(PA), inclination (i), systemic velocity (Uys), and rotation 
velocity (Uot)- These initial estimates were then used as 
inputs to build model-cubes and subsequently refined iter¬ 
atively until we obtained a model-cube that closely matches 
the observed cube. Given the low S/N ratio and possi¬ 
ble contamination from the underlying Hi debris, the best 
models were not determined by a blind minimization, 
but by visual inspection of channel maps and PV-diagrams. 

The intial models are built assuming a thin disc with 
an exponential surface density profile, a flat rotation curve, 
and a constant Hi velocity dispersion (cthi). The disc is 


projected on the sky (according to the values of xq, yo^ PA, 
i, and Vsys) to generate a model-cube, which is subsequently 
smoothed to the same spatial and spectral resolutions of the 
observations. Finally, the model-cube is renormalized on a 
pixel-by-pixel basis to the observed Hi map, i.e. we repro¬ 
duce the asymmetric gas distribution by imposing that, at 
every spatial pixel, the flux along each Hi line profile in 
the model-cube is equal to that of the corresponding line 
profile in the observations (see Lelli et al. 2012a,b, for de¬ 
tails). As a result of this normalization, the assumption of 
an exponential Hi surface density profile has little effect 
on the final model, but is still useful to (i) reproduce beam- 
smearing effects due to a varying HI flux on spatial scales 
smaller than the HI beam, and (ii) obtain an estimate of i 
by comparing model HI maps (before normalization) to the 
observed Hi maps (see dashed ellipses on the top-middle 
panels of Figs. 5 to 10). For VCC 2062, the model-cube was 
not normalized to the observed Hi map because two dis¬ 
tinct kinematic components overlap in both space {R 15" 
to 30" towards the South-West) and velocity (V.o.s. — 1110 
to 1160 kms“^); a simple exponential disc model allows dis¬ 
entangling these two kinematic components (see Fig. 10). 

Once the geometrical parameters of the Hi disc are 
fixed, the free parameters of the model are (Jhi , Vot, and 
the outermost disc radius Rout- These are largely inde¬ 
pendent as they have different effects on the model: Uot 
and Rout increases/decreases the extent of PV-diagrams in 
the vertical and horizontal directions, respectively, whereas 
(Jhi changes the overall shape of PV-diagrams and in¬ 
creases/decreases the distance between density contours. 

5.1.1. Comparing models and observations 

The bottom panels of Figs. 5 to 10 compare PV-diagrams 
from the observed cubes and the model-cubes. For all six 
TDGs, a rotating disc with axisymmetric kinematics can 
reproduce the Hi velocity gradient along the major axis. 
The Hi emission along the minor axis is also well repro¬ 
duced, suggesting that any non-circular motion should be 
relatively small. Some radial motions may be present in 
NGC 5291N and VCC 2062, as hinted to by the tilt of their 
minor-axis PV-diagrams. For these two objects, we built ad¬ 
ditional models including a global radial component (Vad)- 
It is unclear whether the inclusion of Uad substantially im¬ 
proves the models, but we can constrain the ratio Vad/Kot 
to be smaller than 0.3. Hence, the disc kinematics is domi¬ 
nated by rotation and radial motions can be neglected to a 
first approximation. We stress that the inclusion of Uad has 
no effect on the measured values of Vj-ot • The quality of our 
modelling is further demonstrated in Appendix C, where 
we compare HI channel maps from the observed cubes and 
the model-cubes. 

In general, our initial estimates of the kinematical pa¬ 
rameters gave a good description of the observations, but 
we have also built models with different parameter values 
to estimate the corresponding uncertainties. The parame¬ 
ters of the final models are listed in Table 7. For cthi we 
find typical values between 7 and 12 kms“^ apart from 
NGC 5291N and NGC 5291S. For these two TDGs, a ve¬ 
locity dispersion of ^20 kms“^ is required to reproduce the 
observed cube: this high value may be either physical (high 
turbulence in the ISM) or due to unresolved non-circular 
motions within the beam. In NGC 5291 a higher turbulence 
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Fig. 5: Top panels: optical image (left), total FI I map (middle)^ and Hi velocity field {right). The dashed ellipse corre¬ 
sponds to the disc model described in Sect. 5.1. The cross and dashed line illustrate the kinematical centre and major 
axis, respectively. In the bottom-left corner, we show the linear scale (optical image) and the Hi beam (total Hi map 
and velocity field) as given in Table 6. In the total Hi map, contours are at ^4.5, 9, 13.5, 18, and 22.5 Mq pc“^. Bottom 
panels: PV-diagrams obtained from the observed cube (left), model-cube {middle), and residual-cube {right) along the 
major and minor axes. Solid contours range from 2cr to 8cr in steps of la. Dashed contours range from —2a to —4a in 
steps of —la. The horizontal and vertical lines correspond to the systemic velocity and dynamical centre, respectively. 


may be expected because the HI ring is the result of a fast 
head-on collision with a massive companion. 

Although a rotating disc provides a good description 
of the Hi kinematics, there is an important issue with this 
model. The orbital times at the outer edge of the disc (torb) 
are much longer than the estimated ages of the TDGs, im¬ 
plying that the discs would have undergone less than a 


full orbit. In Table 7, we provide both torb and the ra¬ 
tio tmerg/^orb, where tmerg IS the characteristic dynamical 
time of the interaction/merger event. In particular, we as¬ 
sume that (i) the Hi ring around NGC 5291 has formed 
^360 Myr ago (B07), (ii) the progenitor galaxies form¬ 
ing NGC 7252 fell together ^700 Myr ago (Hibbard & 
Mihos 1995; Chien & Barnes 2010), and (hi) the interac- 
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Fig. 6: Same as Fig. 5 for NGC 5291S. In the total Hi map, contours are at ^4.5, 9, 13.5, and 18 Mq pc 
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tion/merger forming VCC 2062 and NGC 4694 has hap¬ 
pened between 0.5 and 1 Gyr ago. Note that tmerg is a hard 
upper limit for the TDG ages, since it takes time for the 
tidal debris to be ejected and for the gas within them to 
collapse under self-gravity. It is unclear, therefore, whether 
the Hi discs associated with the TDGs have had enough 
time to reach dynamical equilibrium. This issue is further 
discussed in Sect. 6.1 with the aid of numerical simulations. 

Motivated by the long orbital times, we have explored 
alternative models and found that a collision between two 
large Hi clouds could also reproduce the Hi data. These 
models, however, need to be highly fine-tuned, given that 
the two Hi clouds must have sizes, spatial separation, and 


velocity separation mimicking the symmetric sides of a ro¬ 
tating disc. This scenario, therefore, appears contrived. We 
conclude that we are likely observing the formation of ro¬ 
tating Hi discs in the nearby Universe. 

5.1.2. Model assumptions and uncertainties 

Dwarf galaxies can show both flat and rising rotation curves 
(e.g. Swaters et al. 2009; Lelli et al. 2012a,b, 2014b; Oh 
et al. 2015). We assume that Uot is constant with R to 
minimize the number of free parameters in our models. This 
assumption has little effects on our results because we are 
interested in the rotation velocity at the outermost radius 
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Fig. 7: Same as Fig. 5 for NGC 5291SW. In the total Hi map, contours are at ^4.5, 9, and 13.5 Mq pc 
Table 7: Properties of the TDGs from 3D kinematical models 


Object 

R.A. 

(J2000) 

Dec. 

(J2000) 

Fsys 

(kms“^) 

PA 

D 

i 

n 

Rout 

(kpc) 

Vrot 

(kms“^) 

cthi 

(kms“^) 

torh 

(Gyr) 

tmerg / Grb 

NGG 5291N 

13 47 20.3 

-30 20 54 

4150 

0 

55 

4.8+1.2 

40+9 

20+2* 

0.7+0.2 

0.5 

NGG 5291S 

13 47 22.7 

-30 27 40 

4660 

5 

65 

7.2+1.2 

20+6 

20+2* 

2.2+0.7 

0.2 

NGG 5291SW 

13 47 17.1 

-30 28 21 

4688 

50 

60 

4.8+1.2 

22+7 

12+2 

1.3+0.4 

0.3 

NGG 7252E 

22 20 56.0 

-24 41 10 

4815 

100 

80 

4.5+1.9 

11+5 

10+2 

2.5+1.6 

0.3 

NGG 7252NW 

22 20 33.7 

-24 37 24 

4678 

125 

60 

7.7+1.9 

16+6 

10+2 

3.0+1.2 

0.2 

VGG 2062 

12 48 00.8 

+10 58 15 

1152 

55 

45 

2.6+0.6 

13+7 

7+1 

1.2+0.5 

0.4-0.8 


Notes. *These high values of cthi may be either physical or due to unresolved non-circular motions within the Hi beam. 
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Fig. 8: Same as Fig. 5 for NGC 7252E. In the total Hi map, contours are at ^1.5, 3.0, and 4.5 Mq pc 


(K)ut)- The actual shape of the rotation curve is difficult 
to determine due to the low number of resolution elements 
(< 4 beams along the disc major axis): in this situation a 
flat rotation curve produces similar model-cubes as a solid- 
body one (cf. Lelli et al. 2012a, Figure 5, second row). As 
a further test, we built additional models assuming a solid- 
body rotation curve (see Appendix D). For NGC 5291N 
and NGC 7252NW, a flat rotation curve is preferable to a 
solid-body one, as indicated by the observed HI emission at 
high rotation velocities near the disc center. For all TDGs, 
solid-body and flat rotation curves return similar values of 
Vout but there is a weak systematic effect leading to higher 
Vout for solid-body rotation curves by r\j 2-3 km s This 


is well within our assumed errors. Hence, the dynamical 
masses estimated in Sect. 5.2 are largely independent of 
the assumed rotation curve shape. 

The uncertainties on (xq, PA, and V^ys have neg¬ 
ligible effect on the measured values of l^ot- PV diagrams 
from model-cubes and observed cubes show that these pa¬ 
rameters have been properly determined (see Figs. 5 to 10, 
bottom panels). The value of i requires special attention 
because Vj-ot oc V[.o.s./sin(i), where Hi.o.s. is the projected 
velocity along the line-of-sight. We constrain the values 
of i by comparing model Hi maps (before pixel-to-pixel 
renormalization) to the observed Hi maps (see dashed el¬ 
lipses in Figs. 5 to 10). In general, we find that inclinations 
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Fig. 9: Same as Fig. 5 for NGC 7252NW. In the total Hi map, contours are at ^1.5, 3.0, 4.5, and 6 Mq pc 
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higher/lower than the assumed ones by ^30° are excluded 
by the data. This is conservatively considered as a 2cr error. 
The error on T^ot is then given by 


^Kot = 


r 


2 


/ 


i 

V + 

/ 

_sin(i) _ 


_ '■°4an(z)_ 


1 2 


( 1 ) 


where = 15° and „ g is estimated as half the velocity 
resolution of the observations. 


5.2. Mass Budget 

In the previous section, we showed that rotating disc mod¬ 
els provide a good description of the HI gas associated with 
TDGs. It is unclear, however, whether these Hi discs are in 
dynamical equilibrium, given that they did not have enough 
time to complete a full orbit since the interaction/merger 
event. In the following, we assume that the TDGs in our 
sample are in dynamical equilibrium and estimate their dy¬ 
namical masses from the Hi rotation velocities. The in¬ 
ferred dynamical masses are then compared with the ob¬ 
served baryonic masses to determine the DM content of 
TDGs. As we discuss in Sect. 6.1, numerical simulations 
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Fig. 10: Same as Fig. 5 for VCC 2062 In the total Hi map, contours are at ^1.5, 3.0, 4.5, 6, and 7.5 Mq pc 
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suggest that these dynamical masses may be quite reliable 
despite the long orbital times. 

5.2.1. Dynamical mass 

For the TDGs in our sample, the rotation velocities (Vot) 
are comparable with the Hi velocity dispersions (cthi), 
thus the former must be corrected for pressure support in 
order to trace the gravitational potential. We apply the 
asymmetric-drift correction using the following equation: 

Virc = V^3t + 4l(^AHl), (2) 

where Vdrc is the circular velocity tracing the potential well. 
This equation assumes that (i) the velocity dispersion is 


isotropic, (ii) the velocity dispersion and scale height of 
the disc are constant with radius, and (iii) the Hi surface 
density profile is an exponential function with scale length 
hui (see Appendix A of Lelli et al. 2014a, for details). 
Using 3D disc models, we estimate that our TDGs have 
Rout/hui — 2, where Rout is the outermost radius of the 
Hi disc. The dashed ellipses in Figs. 5 to 10 (top-middle 
panels) show the model Hi maps at R 2hHi (before 
being renormalized to the observed HI map, see Sect. 5.1). 
The only exception is NGC 5291N, which have Rout — ^hi • 

The dynamical mass within Rout is then given by 

Mdyn = e (3) 
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Table 8: Mass budget in TDGs assuming dynamical equilibrium 


Object 

Fcirc 

(kms“^) 

Mdyn 
(10® Mo) 

Afatom 
(10® Mo) 

M. 

(10® Mo) 

Mmol 
(10® Mo) 

Mbar 

(10® Mo) 

Mdyn/Mbar 

NGG 5291N 

45±9 

23±11 

12.3±1.2 

1.1±0.5 

2.2±1.1 

15.6±1.7 

1.5±0.7 

NGG 5291S 

35±6 

20±8 

11.8±1.2 

0.8±0.4 

3.3±1.7 

15.9±2.1 

1.3±0.5 

NGG 5291SW 

28±7 

8.7±4.9 

4.9±0.5 

0.3±0.2 

1.9±1.0* 

7.1±1.2 

1.2±0.7 

NGG 7252E 

18±5 

3.4±2.4 

2.9±0.3 

0.85±0.43 

0.20±0.10* 

3.9±0.5 

0.9±0.6 

NGG 7252NW 

21±6 

7.9±4.9 

7.0±0.7 

0.94±0.47 

0.22±0.11 

8.2±0.8 

1.0±0.6 

VGG 2062 

16±7 

1.5±1.3 

0.75±0.07 

0.47±0.24 

0.23±0.12 

1.4±0.3 

1.0±0.9 


Notes. *For these two galaxies GO observations are not available, thus Mmoi is indirectly estimated from the SFR assuming that 
TDGs have the same SFE as spiral galaxies 5 x 10“^° yr~^). See Sect. 5.2.2 for details. 


where G is Newton’s constant and 5 is a parameter of the 
order of unity that depends on the 3D distribution of the 
total mass. For simplicity, we assume e = 1 (spherical ge¬ 
ometry). For a thin exponential disc, one can calculate that 
5 0.95 8 it R = h and 5 0.77 ai R = 2h (cf. Freeman 

1970). The resulting dynamical masses are given in Table 8. 
Besides the uncertainty due to the dynamical equilibrium 
of TDGs, the uncertainty on M^yn is driven by errors on 
Vcirc and i?out- The former error is estimated using Eq. 1, 
while the latter is conservatively estimated as half the spa¬ 
tial resolution of the HI data. 

5.2.2. Baryonic mass 

The baryonic mass of a galaxy (Mbar) is given by three main 
components (e.g. McGaugh 2012): atomic gas (Matom), 
molecular gas (Mmol), and stars (M*). In TDGs Matom 
largely dominates the baryonic mass budget, thus the un¬ 
certainties on M* and Mmol due to the choice of, re¬ 
spectively, a stellar mass-to-light ratio (T*) and CO-to- 
H2 conversion factor (Xco) have minimal impact. On the 
other hand, measuring Matom is not straightforward be¬ 
cause TDGs are still embedded within gas-rich tidal debris. 
It is reasonable to assume that all the HI gas within Rout is 
associated with the TDGs, since they represent a local po¬ 
tential well that “clears out” the gas from the surroundings. 
The values of Matom in Table 8 are based on this assump¬ 
tion and adopt a typical error of 10% (the Hi masses are 
multiplied by a factor of 1.37 to take the contributions of 
Helium and heavier elements into account). Alternatively, 
one may assume that, at the TDG location, some gas is 
associated with the underlying tidal debris. In the extreme 
hypothesis that tidal debris has a uniform HI surface den¬ 
sity of 1 Mq pc“^ over the TDG area, the values of Matom 
in Table 8 would decrease by ^30% to 40%. 

As we discussed in previous sections, TDGs have higher 
metallicities than typical dwarfs, thus CO lines are readily 
detected and H 2 masses can be estimated using a typical 
value of Xco = 2.0 x 10^^ cm^ (K kms“^)“^. For three 
TDGs (NGC 5291N, NGC 5291S, and NGC 7252NW), we 
use the molecular masses from Braine et al. (2001), rescaled 
to the distances in Table 1. For VCC 2062, we adopt the 
molecular mass from D07. For the remaining two TDGs 
(NGC 5291SW and NGC 7252E), CO observations are not 
available, thus we estimated Mmoi assuming that TDGs 
and spiral galaxies have the same star-formation efficiency: 
SFE = SFR/Mmoi — 5 X 10“^^ yr“^ (Leroy et al. 2008). 
This is observationally motivated by the results of Braine 


et al. (2001) for a sample of eight TDGs. For these two 
TDGs, SFRs are estimated using FUV fluxes from Boquien 
et al. (2009) and the calibration from Kennicutt & Evans 
(2012), which is based on Starburst99 models (Leitherer 
et al. 1999) and a Kroupa initial mass function (Kroupa 
& Weidner 2003). For all TDGs, we adopt a conservative 
error of 50% on Mmol to consider possible variations in Xqo 
and/or SFE among galaxies. The molecular masses of these 
TDGs are from ^4 to ^30 times smaller than Matom, thus 
they marginally contribute to Mbar- 

The stellar masses of TDGs are tricky to estimate be¬ 
cause they are formed by a complex mix of young stars 
(comprising the dominant stellar population) and possibly 
older stars ejected from the parent galaxies. For the three 
TDGs around NGC 5291, we adopt the stellar masses from 
B07, which were derived using iF-band luminosities and 
adopting T* = 0.3. For the other three TDGs, we adopt 
the values from Table 4 of Due et al. (2014), which were 
obtained by re-analysing the photometric data of Boquien 
et al. (2009, 2010) and fitting spectral-energy distributions 
with stellar population models. We adopt a conservative 
error of 50% on M^. The stellar masses of these TDGs are 
from ^3 to ^20 times smaller than Matom, thus have little 
impact on Mbar- The only exception is VCC 2062 that has 
Matom — 1-5M^ and represents the oldest and most evolved 
TDG in our sample. 

The mass budget of our six TDGs is summarized in 
Table 8. Assuming that these TDGs are in dynamical equi¬ 
librium, we find that their dynamical masses are fully con¬ 
sistent with the observed baryonic masses, implying that 
they are nearly devoid of DM. Thus, we do not confirm the 
results of B07, who reported values of M^yn/Mbar between 
2 and 3 for three TDGs around NGC 5291. The reasons 
for this discrepancy are discussed in Appendix B. We note 
that the ratio M^yn/Mbar depends on the assumed distance 
as 1 /D. NGC 4068 and VCC 2062 lie in the Virgo Cluster: 
assuming D = 17 ± 1 Mpc (cf. Mei et al. 2007), the result¬ 
ing error on M^yn/Mbar is smaller than 6%. NGC 5291 and 
NGC 7252 are relatively far away (Vys ^ 4300 km s“^) and 
the Hubble flow is a reliable distance indicator: even assum¬ 
ing an extreme peculiar velocity of 500 km s“^ along the 
line of sight, the resulting error on Mdyn/Mbar is smaller 
than 12%. 

5.3. Comparison with other dwarf galaxies 

To assess whether the DM content of TDGs is really anoma¬ 
lous, we use a control-sample of typical star-forming dwarfs 
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Fig. 11: Mass budget within 3 h* in different types 

of dwarf galaxies: Irrs (blue diamonds, from Swaters et al. 
2009), BCDs (black squares, from Lelli et al. 2014b), and 
TDGs (red points). For Irrs Mbar is estimated assuming 
both a “heavy disc” with = 2 (open symbols) and a 
“light disc” with Tf = 1 (filled symbols). For BCDs the 
stellar mass is estimated from color-magnitude diagrams of 
resolved stellar populations (see Lelli et al. 2014b). 


formed by 30 irregulars (Irrs) and 6 blue compact dwarfs 
(BCDs). The rotation velocities of these galaxies have been 
measured by Swaters et al. (2009) and Lelli et al. (2014b) 
building model-cubes of the Hi observations, thus they are 
fully comparable with those derived here. TDGs are heavily 
gas-dominated (Mgas/M* 10), thus Mdyn/Mbar is esti¬ 
mated at roughly R 2hHi, where hni is the scale-length 
of the Hi disc (Sect. 5.2). To have a meaningful compari¬ 
son with typical dwarfs, we need to compute Mdyn/Tfbar 
within a similar radius. In general, these galaxies have 
Mgas/M^ <1 — 2, thus it is more natural to consider the 
scale-length of the stellar disc h*. Given that typical dwarfs 
have hni/h* 1.5 (cf. Swaters et al. 2002), we estimate 
Mdyn/Mbar 8 it R 3h^ (the edge of the stellar disc). 

The Irrs are drawn from Swaters et al. (2009) using the 
following criteria: (i) the galaxies have 30° ^ i ^ 80°, (ii) 
the Hi rotation curves have quality flag q > 2 and are 
traced out to R 3h*, and (hi) Vrot < 100 kms“^ at the 
last measured point. The BCDs are drawn from Lelli et al. 
(2014b) considering only galaxies with reliable estimates of 
Vrot (see their Sect. 7). 

Contrary to TDGs, typical dwarfs often have com¬ 
parable stellar and gaseous masses, thus the values of 
Mdyn/Mbar strougly depend on the assumed stellar mass- 
to-light ratio T>^. For Irrs we compute two different values 
of Mdyn/Tfbar assumiug a “heavy disc” (Tf = 2) and a 
“light disc” (Tf = 1). We do not provide an error bar 
on Mdyn/Tfbar bccausc this is fully dominated by the as¬ 
sumed variation in T^. For BCDs, we use the stellar masses 
obtained from modelling color-magnitude diagrams of re¬ 
solved stellar populations (see Lelli et al. 2014b, for de¬ 


tails and references). Our results are illustrated in Fig. 11. 
Despite the large error bars, it is clear that TDGs sys¬ 
tematically have smaller values of Mdyn/Mbar than typical 
dwarfs, pointing to an apparent lack of DM. 


6. Discussion 

In this paper we studied the gas dynamics in a sample of six 
TDGs. These galaxies have higher metallicities than dwarfs 
of similar mass and appear kinematically decoupled from 
the surrounding tidal debris, indicating that they are bona- 
fide TDGs (Sect. 4). We built 3D kinematical models and 
found that the Hi emission associated with TDGs can be 
described by rotating discs (Sect. 5.1). The orbital times 
of these discs, however, are much longer than the merger 
timescales (as inferred from numerical simulations), indi¬ 
cating that they have undergone less than a full orbit since 
the TDG formation. This raises the question as to whether 
TDGs are in dynamical equilibrium and the observed ve¬ 
locities can be used to estimate dynamical masses. If one 
assumes that these TDGs are in dynamical equilibrium, 
the inferred dynamical masses are fully consistent with the 
observed baryonic masses, implying that TDGs are nearly 
devoid of DM (Sect. 5.2). In the following, we start by scru¬ 
tinizing the hypothesis of dynamical equilibrium with the 
aid of numerical simulations. Then, we discuss the impli¬ 
cations of our results for the distribution of non-baryonic 
DM in galaxies and the possible amount of dark baryons. 
Finally, we investigate the location of TDGs on the bary¬ 
onic Tully-Fisher relation (BTFR, McGaugh et al. 2000) 
and discuss alternative theories like MOND. 

6.1. Dynamical equilibrium: clues from numerical simulations 

Bournaud & Due (2006) and B07 presented hydro- 
dynamical simulations of interacting/merging systems, in¬ 
vestigating the detailed process of TDG formation. We have 
re-analysed the output of these simulations to check the hy¬ 
pothesis that the relaxation time in TDGs could be shorter 
than the orbital time. For different snapshot of the simu¬ 
lations, we measured (i) the dynamical mass M^yn of the 
simulated TDGs inferred from gas rotation curves, i.e. con¬ 
sidering the rotation velocities of individual gas particles 
with respect to the mean spin axis, and (ii) the real mass 
Mreai within the TDG radius, obtained by summing gas, 
star, and DM particles (the latter ones constitute a negligi¬ 
ble amount of mass). For a galaxy in dynamical equilibrium, 
we expect Mdyn/Mreai = 1 within ~10% (the typical errors 
associated with the estimate of M^yn from simulated gas ro¬ 
tation curves). In Fig. 12, the ratio M^yn/Mreai is plotted 
versus the look-back time, where t = 0 corresponds to the 
snapshot of the simulation that reproduces the present-day 
properties of the observed system. We show two different 
simulations: red dots are TDGs formed in the tidal tails of 
a standard disc-disc merger (similar to NGC 7252), black 
crosses are TDGs formed in a collisional ring (matching 
the properties of NGC 5291). These simulated TDGs have 
dynamical masses comparable to those of observed TDGs 
(cf. B07). At early times, the values of M^yn/Mreai signifi¬ 
cantly differ from 1 and display a large scatter, indicating 
that TDGs are still far from dynamical equilibrium. After 
a few 100 Myr, however, the scatter substantially decrease 
and the values of Mdyn/Mreai gets closer and closer to 1. 
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the Milky Way and found no evidence for an accreted stel¬ 
lar disc, implying that near-plane mergers have been very 
rare in our Galaxy during the past ^10 Gyr. Nonetheless, 
a “dark disk” may help to explain the vertical kinematics 
of stars in the solar neighborhood (Bienayme et al. 2014). 
Moreover, Fan et al. (2013) and McCullough & Randall 
(2013) speculated that exotic forms of self-interacting DM 
can form rotating dark discs as massive as the baryonic 
ones. The values of M^yn/Mbar in TDGs shed new light on 
these putative dark discs, given that any kind of rotation- 
supported material in the progenitor galaxies may form 
tidal tails and be accreted in TDGs. Our results imply 
that either these putative dark discs are significantly less 
-250 -200 -150 -100 -50 0 massive than the observed baryonic discs or they are too 

Look-back time (Myr) dynamically-hot to be accreted by TDGs. 


I .U 
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Fig. 12: The ratio Mdyn/Tfreai versus look-back time for 
simulated TDGs. Mdyn is the dynamical mass inferred from 
the mean rotation velocity of gas particles, while Mreai is 
the mass obtained by adding gas, stars, and DM particles 
within the TDG radius. The time t = 0 corresponds to the 
snapshot of the simulation that reproduces the present-day 
large-scale properties of the interacting system. Red dots 
are TDGs formed in the tidal tails of a standard disc-disc 
merger that resembles NGC 7252. Black crosses are TDGs 
formed in a collisional HI ring that matches the properties 
of NGC 5291 (B07). 

This suggests that the simulated TDGs may have reached 
dynamical equilibrium in less than an orbital time. 

For collisionless stellar systems, it is known that the 
process of relaxation takes a few orbital times, since it oc¬ 
curs through the energy exchange between individual par¬ 
ticles and the total gravitational potential evolving over 
orbital timescales (Lynden-Bell 1967). For collisional gas- 
dominated systems (such as TDGs), the energy can be ex¬ 
changed and radiated away on a local crossing time of the 
order of ^10 Myr (Mac Low 1999). Thus, one may speculate 
that in collisional gas-rich systems the process of relaxation 
takes less than one orbital time and is significantly faster 
than in collisionless ones. Further theoretical studies are 
needed to understand these processes. It is clear, however, 
that TDGs provide a unique opportunity to study galaxy 
relaxation and represent a local benchmark for understand¬ 
ing disc formation in gas-rich galaxies at high redshifts. 

6.2. Implications for the DM distribution 

Assuming that the TDGs in our sample are in dynamical 
equilibrium, we deduce that they are devoid of DM as ex¬ 
pected from numerical simulations where tidal forces segre¬ 
gate disc material from halo material (Barnes & Hernquist 
1992; Bournaud & Due 2006; Wetzstein et al. 2007). This 
provides an indirect but empirical confirmation that most 
DM must be distributed in pressure-supported halos. 

Using cosmological simulations. Read et al. (2008, 2009) 
and Pillepich et al. (2014) argued that spiral galaxies like 
the Milky Way should possess a rotating disc of non- 
baryonic DM. This forms because merging satellites are 
preferentially dragged towards the disc plane by dynami¬ 
cal friction and then disrupted by tides, depositing both 
stars and DM into the disc. Ruchti et al. (2014, 2015) used 
a chemo-dynamical template to search for accreted stars in 


6.3. Implications for dark gas 

Several observations suggest that a fraction of gas in galaxy 
discs is not accounted for by Hi or CO observations. In 
the Milky Way, this “dark gas” has been inferred from (i) 
7 -ray emission due to the collision of cosmic rays with in¬ 
terstellar nucleons (Grenier et al. 2005; Abdo et al. 2010), 
and (ii) a measured excess of thermal dust emission with 
respect to a canonical dust-to-gas ratio, as indicated by 
FIR/submm maps from IRAS (Reach et al. 1994), COBE 
(Reach et al. 1998), and Planck (Planck Collaboration et al. 
2011). Both techniques indicate that, in the Milky Way, the 
dark gas mass is comparable to the CO-traced molecular 
mass. The relative mass of dark gas may be even larger in 
the LMC (Galliano et al. 2011). This dark gas may con¬ 
sist of optically thick Hi and/or H 2 in the outer regions 
of molecular clouds, where CO is photodissociated (Wolfire 
et al. 2010; Glover et al. 2010). TDGs offer a unique way to 
probe dark gas in galaxies outside the Local Group, given 
that this undetected baryonic component would give values 
of Mdyn/A/bar > 1- Unfortunately, our current estimates of 
Mdyn/Mbar are not precise enough to identify any potential 
contribution from dark gas at the expected level. 

Pfenniger et al. (1994) and Pfenniger & Combes (1994) 
proposed very cold molecular gas as an alternative to non- 
baryonic DM in galaxies. This putative dark gas would 
consist of small “clumpuscules” that do not emit any de¬ 
tectable line and are not associated with star formation. 
Intriguingly, the observed rotation curves of disc galaxies 
can be well reproduced simply by scaling up the observed 
contributions of stars and gas (Hoekstra et al. 2001; Swaters 
et al. 2012). This very cold molecular gas has also been 
supported by other dynamical arguments such as the exis¬ 
tence of spiral patterns and long-lived warps in the outer 
regions of Hi discs (Bureau et al. 1999; Masset & Bureau 
2003; Revaz & Pfenniger 2004; Revaz et al. 2009). These H 2 
clumpuscules should end up in TDGs and be detected as a 
significant excess of dynamical mass over the baryonic one. 
In contrast with B07, our current results do not support 
the existence of large numbers of cold H 2 clumpuscules in 
excess of the more modest expectation for “standard” dark 
gas. 

6.4. The baryonic Tully-Fisher relation and MONO 

Fig. 13 shows the location of TDGs on the BTFR. We use 
data from McGaugh (2005) for star-dominated galaxies and 
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Table 9: MOND analysis assuming dynamical equilibrium. 


Object 

(kpc) 

gNe/do 

gNi/ao 

l^ISO(l) 

(kms-1) 

EiSO(2) 

(kms“^) 

Eefe(i) 

(kms“^) 

VeFE(2) 

(kms-b 

kcirc 

(kms-b 

NGG 5291N 

65 

0.021 

0.073 

77±2 

73±2 

62±7 

57±7 

45±9 

NGG 5291S 

60 

0.025 

0.033 

76±3 

73±2 

54±8 

49±8 

35±6 

NGG 5291SW 

75 

0.016 

0.033 

62±3 

59±3 

46±7 

43±7 

28±7 

NGG 7252E 

50 

0.044 

0.019 

52±2 

50±2 

30±5 

28±5 

18±5 

NGG 7252NW 

80 

0.017 

0.015 

63±2 

61±2 

41±6 

39±6 

21±6 

VGG 2062 

20 

0.024 

0.022 

41±2 

39±2 

27±5 

25±5 

16±7 


Notes. Dp is the projected distance of the TDG to the central galaxy remnant, which is used to estimate ^Ne- 



Fig. 13: The location of TDGs on the BTFR. Vdrc is the cir¬ 
cular velocity after correcting for pressure support (Eq. 2). 
For star-dominated galaxies (McGaugh 2005) and gas- 
dominated ones (McGaugh 2012), Ecirc is measured along 
the flat part of the rotation curve. For TDGs, the shape of 
the rotation curve is uncertain, thus Vdrc iiiay not corre¬ 
spond to the velocity along its flat part. 


McGaugh (2012) for gas-dominated ones. TDGs are sys¬ 
tematically shifted to the left side of the BTFR. This is 
not fully surprising since these TDGs seem to have smaller 
DM content than typical dwarf galaxies (Fig. 11), hence 
they have lower rotation velocities for a given Mbar- Note 
that McGaugh (2005, 2012) used the asymptotic velocity 
along the flat part of the rotation curve (Rflat), which is 
known to minimize the scatter on the BTFR (Verheijen 
2001; Noordermeer & Verheijen 2007). For the TDGs in 
our sample, the shape of the rotation curve is uncertain, 
thus it remains unclear whether we are really tracing Viat- 
The Hi data, however, are consistent with a flat rotation 
curve (Sect. 5.1). 


The BTFR is one of the major successful predictions of 
MOND (e.g. McGaugh 2005, 2012). According to MOND, 
disc galaxies can deviate from the BTFR in only two cir¬ 
cumstances: (i) if they are out of dynamical equilibrium 
(e.g. McGaugh & Wolf 2010) and/or (ii) if they are affected 
by the external-field effect (EFE, Bekenstein & Milgrom 
1984). The numerical simulations discussed in Sect. 6.1 sug¬ 
gest that TDGs may reach dynamical equilibrium after only 
a few 100 Myr (see Eig. 12), but it remains unclear whether 
this would still hold in MOND. Hydrodynamical simula¬ 
tions in MOND are needed to clarify this issue. 

If MOND is interpreted as modified gravity, the strong 
equivalence principle is violated and the internal dynam¬ 
ics of a system can be affected by an external gravita¬ 
tional field (Bekenstein & Milgrom 1984). We define ^Ni = 
GMbar/7?out internal Newtonian acceleration of the 

TDG. We also define ^Ne = GMhost/^p ^.s the external 
Newtonian acceleration, where Mhost is the baryonic mass 
of the central galaxy remnant and Dp is the projected dis¬ 
tance to the TDGs. Mhost is estimated as Lk x T^, where 
I/K is given in Table 1 and = 0.6 Mq/Lq (McGaugh 
& Schombert 2015). The values of ^Ni and ^Ne are given in 
Table 9. They are roughly comparable, hence the EEE may 
be important. 

To a first approximation, the EEE can be described by 
Eq. 60 of Eamaey & McGaugh (2012)^: 


di 


, ^Ni + ^'Ne . , 

gml^[ - ) +^Ne 

do 


. gNi + ^Ne\ 

iy[ - -ly 

ao ) 


gNe \ 

do J \ ’ 

(4) 


where Ui = V^/R is the internal acceleration of the TDG, 
ao — 1.3 X 10“^ cm s“^ is the MOND acceleration scale, 
and u{y) is the interpolation function. Eollowing McGaugh 
(2008), we assume i^niy) = 1/[1—exp(—where the 
cases n = 1 and n = 2 nearly correspond to the so-called 
“simple” and “standard” functions, respectively. In Table 9 
we provide the velocities expected in MOND for n = 1 
and n = 2, considering both the isolated case (Viso with 
^Ne = 0) and the EEE (Vefe)- Eormal errors are estimated 
considering the uncertainties on Mbar and allowing Mhost to 
vary by a factor of 2. The MOND velocities are decreased by 
the EEE, but they are still systematically higher than the 

^ This equation derives from the quasi-linear formulation of 
MOND (Milgrom 2010). We have also considered the non-linear 
formulation of Bekenstein V Milgrom (1984) using Eq. 59 of 
Famaey V McGaugh (2012) and adopting the simple interpola¬ 
tion function. We find consistent results within ~1 kms“^. 
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observed ones. For individual TDGs, the tension between 
MOND and observations is between and ^2cr. 

We stress that Eq. 4 grossly simplifies the EFE because 
it formally applies to the one-dimensional case and neglects 
the direction of the external field. We are also neglecting the 
time evolution of the EEE, which may be relevant here be¬ 
cause TDGs start to form closer to the host galaxy (where 
the EEE is stronger) and are subsequently ejected to larger 
distances at high speeds. Since TDGs have long orbital 
timescales, they may respond in a non-adiabatic fashion 
to the rapidly varying EEE, keeping memory of the ini¬ 
tial, stronger EEE. It would be interesting to investigate 
whether a more rigorous treatment of the EEE (using nu¬ 
merical simulations) could further decrease the velocities 
expected in MOND. 

Einally, we note that MOND can also be interpreted 
as modified inertia (Milgrom 1994, 2006). In this case, the 
theory becomes strongly non-local and the force acting on 
a test particle depends on the full orbital history. A simple 
solution can be obtained in the case of stable, stationary, 
circular orbits (Milgrom 1994), but it is unclear what to ex¬ 
pect for discs that have undergone less than one rotation. 
In this respect, TDGs greatly differ from typical disc galax¬ 
ies and could, in principle, provide a test for such non-local 
theories. 

7. Conclusions 

We presented a systematic study of the gas dynamics in 
a sample of six TDGs around three interacting/merging 
systems: NGC 4694, NGC 5291, and NGC 7252 (“Atoms 
for Peace”). Eor NGC 4694 and NGC 5291, we re-analysed 
existing HI datacubes from D07 and B07, respectively. Eor 
NGC 7252, we presented new Hi observations from the 
Jansky VLA, lEU observations from VLT-GIRAEEE, and 
long-slit spectroscopy from VLT-EORSl. Our results can 
be summarized as follows: 

1. The putative TDGs near NGC 7252 have near solar 
metallicities (12+log(0/H)^8.6), indicating that they 
are currently forming out of pre-enriched tidal material. 
Moreover, they are associated with steep Hi velocity 
gradients that represent kinematic discontinuities along 
the tails, pointing to local potential wells. Similar re¬ 
sults have been found for three TDG candidates near 
NGC 5291 (B07) and one near NGC 4694 (D07). This 
confirms that the six galaxies in our sample are bona- 
fide TDGs. 

2. The Ha kinematics of NGC 7252NW is complex and 
characterized by very broad line profiles. Most likely, 
we are observing H ll regions dominated by turbulence 
and/or non-circular motions due to stellar feedback. 

3. The Hi emission associated with TDGs can be de¬ 
scribed by 3D disc models. These discs, however, have 
orbital times ranging from ~1 to ^3 Gyr, which are 
significantly longer than the TDG formation timescales 
(<1 Gyr). This raises the question as to whether TDGs 
have had enough time to reach dynamical equilibrium. 
TDGs appear to be rotation-supported galaxies that are 
forming at z 0 and, thus, are key objects to study the 
processes of disc building and relaxation. 

4. If one assumes that TDGs are in dynamical equilibrium 
(as suggested by numerical simulations), the inferred 
dynamical masses are fully consistent with the observed 


baryonic masses. This implies that TDGs are devoid 
of DM, as predicted by theoretical models where tidal 
forces segregate baryons in the disc from DM in the halo. 
This also puts constraints on putative dark discs (either 
baryonic or non-baryonic) in the progenitor galaxies. 
Moreover, TDGs seem to systematically deviate from 
the baryonic Tully-Eisher relation. These results provide 
a challenging test for theories like MOND. 
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Appendix A: Comparison with archival VLA data 
and single-dish observations for NGC 7252 

For a sanity check, we retrieved from the VLA archive 
the Hi observations of NGC 7252 from Hibbard et al. 
(1994) and compare them with the new ones. We reduced 
the archival observations following standard procedures in 
AIPS and created a datacube with spatial resolution of 25" 
(robust weights) and velocity resolution of 10.6 km s“^. The 
D-array data were interpolated to match the spectral reso¬ 
lution of the C-array data. 



Fig. A.l: PV-diagram for NGC 7252 obtained by integrat¬ 
ing the Hi emission along the declination axis. The new 
data are shown in grayscale. Contours are at 18, 36, and 72 
mJy/beam. The archival data (published by Hibbard et al. 
1994) are represented with a red contour at 18 mJy/beam. 
Both datasets were spatially smoothed to 25" resolution. 

Fig. A.l shows a PV-diagram obtained by collapsing 
the final datacube along the declination axis, giving a com¬ 
plete overview of the kinematic structure of the system. The 
comparison between the new (grey scale) and archival (red 
contour) data shows some differences on small spatial scales 
due to different noise levels in the two data sets. Overall the 
new VLA data are consistent with the archival ones. 

From the new VLA data we derive a total Hi flux 
of 5.0 Jy kms“^. This is in good agreement with the 
single-dish measurement by Richter et al. (1994), who find 
a total Hi flux of 4.6 Jy kms“^ using the Green Bank 
Telescope (FWHM = 21' x 21'). Other single-dish mea¬ 
surements (who rely on single pointings) reported lower 
fluxes: Huchtmeier (1997) find 3.8 Jy kms“^ using the 
Effelsberg telescope (FWHM = 9' x 9'), while Dupraz 
et al. (1990) find 3.6 Jy kms“^ using the Nancay telescope 
(FWHM = 4' [E - W] X 22' [N - S]). Most likely, the lat¬ 
ter single-dish observations miss part of the HI emission in 
the tidal tails, extending for ^12' in the E-W direction. We 
conclude that no significant emission is missing from the 
new VLA data. 


Appendix B: Comparison with B07 for NGC 5291 

Eor the TDGs around NGC 5291, we find significantly lower 
values of M^yn/ATbar than B07 because we estimate both 
higher values of Mgas (by factors of ^1.5 to ^2) and lower 
values of Vot (by a factor of ~2). There are two reasons 
for the different estimates of Mgas^ (i) we integrate the Hi 
flux out to larger radii than B07, according to the results of 
our 3D disc models; this is the dominant effect; and (ii) we 
consider all the Hi emission that is consistent with a ro¬ 
tating disc, whereas B07 assumed that some Hi gas at the 
high/low velocity end was in the background/foreground 
(instead of being associated with the TDG potential well). 
The different estimates of Vot are also driven by two sepa¬ 
rate reasons: (i) we find that line-of-sight velocities at the 
edges of the PV-diagrams (employed by B07) tend to over¬ 
estimate the “true” rotation velocity because they do not 
consider the line-broadening due to turbulent motions; and 
(ii) we find higher disc inclinations than assumed by B07. 
Hereafter, we discuss these two effects in details. 

B07 derived rotation curves using the envelope-tracing 
method, which considers line-of-sight velocities near the 
edges of PV-diagrams. Specifically, for a given radius along 
the major axis, they took the velocity towards the far edge 
of the PV-diagram at 50% of the HI peak level. Eor poorly- 
resolved rotating discs, this is a rough way to consider 
beam-smearing effects, which are well-known to artificially 
broaden the emission-line profiles and systematically skew 
them towards the systemic velocity. The details of beam- 
smearing effects, however, depend on the intrinsic gas dis¬ 
tribution, rotation curve shape, gas velocity dispersion, and 
disc inclination. Thus, full 3D models are required to take 
them into account (Swaters et al. 2009; Lelli et al. 2010). 
In particular, the velocity at 50% of the Hi peak may sys¬ 
tematically over-estimate the rotation velocity in discs with 
low values of Kot/<^Hi (as dwarf galaxies), given that the 
broadening of the line profile is driven by both cthi and 
resolution effects. Our 3D disc models take all these obser¬ 
vational effects into account and suggest that the rotation 
velocities of these TDGs are closer to the Hi peak than 
adopted by B07. This is demonstrated in Eig. B.l. 

Eor the three TDGs near NGC 5291, we built addi¬ 
tional 3D models following the same procedures of Sect. 5.1 
and adopting the kinematical parameters of B07 {i = 45°, 
(Jhi =10 km s“^, and Vot from 50 to 70 km s“^ depending 
on the TDG). The right panels of Eig. B.l compare PV- 
diagrams along the major axis obtained from the observed 
cubes and these model-cubes. The kinematic parameters 
from B07 produce PV-diagrams that are too extended in 
the velocity direction. This leads to large negative residuals 
at high/low line-of-sight velocities and large positive resid¬ 
uals near the systemic velocity, indicating that Vot sin(i) is 
overestimated and am is underestimated. It is clear that 
our new kinematic parameters provide a much better de¬ 
scription of the data (cf. with Eigs. 5, 6, and 7). 

Regardless of the method employed to estimate the ro¬ 
tation velocities, these projected quantities must be cor¬ 
rected for inclination by multiplying by 1/sin(i). B07 as¬ 
sumed that the TDGs around NGC 5291 have the same 
inclination angle as the collisional Hi ring {i = 45°), which 
was estimated using numerical simulations. B07 pointed out 
that this value of i is uncertain because (i) the modelled 
ring is not exactly circular and its morphology depends on 
the details of the collision, and (ii) some simulated TDGs 
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Fig. B.l: Disc models for NGC 5291N (top), NGC 5291S {middle)^ and NGC 5291SW (bottom) adopting the parameters 
from B07. Left panels: total Hi maps from the observed cube (solid contours) and the model-cube (dashed ellipse). The 
cross and dashed line illustrate the kinematical centre and major axis, respectively. The circle to the bottom-left corner 
shows the Hi beam. Right panels: PV diagrams along the major axis obtained from the observed cube, model-cube, and 
residual-cube. Solid contours range from 2cr to 8cr in steps of la. Dashed contours range from —2a to —4a in steps of 
— la. The horizontal and vertical lines correspond to the systemic velocity and dynamical centre, respectively. 


show a discrepancy between their rotation axis and the ring 
axis (up to 18° in a few cases). Using 3D disc models, we 
directly estimate i by comparing observed Hi maps with 
model Hi maps (before pixel-to-pixel renormalization). In 
left panels of Fig. B.l, the dashed ellipses correspond to disc 
models with i = 45°. This value of i is roughly acceptable 
for NGC 5291N and NGC 5291SW, whereas it is clearly 
underestimated for NGC 5291S. The dashed ellipses in the 
top-middle panels of Figs. 5, 6, and 7 illustrate our adopted 
inclination, which are slightly larger than 45° (from 55° to 
65°, as given in Table 7). Thus our estimates of Uot are 
further decreased with respect to those of B07 by a fac¬ 
tor sm(i = 45°)/sin(inew): be. by ^15% to ^30%. Clearly, 
inclination plays a minor role here. 

Appendix C: Channel Maps of individual TDGs 

Figures C.l and C.2 show Hi channel maps from the ob¬ 
served cubes (red contours) and model-cubes (blue con¬ 


tours). Contours are at -3a (dashed), 3a, 6a, and 9a. The 
channel maps are superimposed on an optical image of the 
TDG. The cross marks the kinematical center. Line-of-sight 
velocities are indicated to the top-left corner. The HI beam 
is shown to the bottom-left corner. The model-cubes are 
described in detail in Sect. 5.1. Here we provide a concise 
description for each TDG. 

NGC 5291N is surrounded by Hi emission belonging to the 
underlying tidal debris. At U.o.s. — 4145 to 4100 kms“^, 
there is extended Hi emission to the West, which may cor¬ 
respond to another, more uncertain TDG candidate. 

NGC 5291S is nicely reproduced by our disc model. At ap¬ 
proaching velocities (U.o.s. — 4645 — 4635 kms“^), there is 
some anomalous gas on the receding side of the disc. 

NGC 5291SW is a complex case due to the low signal-to- 
noise ratio. The Hi emission at U.o.s. — 4730 — 4720 kms“^ 
does not show a continuous, coherent kinematical structure, 
hence it has not been considered in the disc model. 

NGC 7252E is well reproduced by our disc model despite 
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it is a poorly resolved case with low signal-to-noise ratio. 
NGC 7252NW is closely reproduced by our disc model. 
VCC 2062 is characterized by two distinct kinematic com¬ 
ponents, which overlap in both space and velocity (see PV- 
diagrams in Fig. 10). The second, irregular component ex¬ 
tends to the South-West at fb.o.s. — 1155 — 1125 kms“^ and 
does not show a coherent velocity structure, hence it is not 
considered in our disc model. 

Appendix D: Additional 3D disc models 

Figures D.l and D.2 compare 3D disc models based on a 
flat rotation curve (middle panel) and a solid-body rotation 
curve (right panel). The models are described in detail in 
Sect. 5.1. The PV-diagrams are obtained along the disc 
major axis. Solid contours range from 2cr to 8cr in steps of 
Icr. Dashed contours range from —2cr to — 4cr in steps of 
— Icr. The horizontal and vertical lines correspond to the 
systemic velocity and dynamical centre, respectively. 
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Fig. C.l: Hi channel maps for NGC 5291N (top), NGC 5291S (middle), and NGC 5291SW (bottom). See Appendix C 
for a detailed description of this image. 
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Fig. C.2: Hi channel maps for NGC 7252E (middle), and NGC 7252NW (bottom), and VCC 2062 (top). See Appendix 
C for a detailed description of this image. 
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Fig. D.l: Major axis PV-diagrams for NGC 5291N (top), NGC 5291S (middle), and NGC 5291SW (bottom). See Appendix 
D for a detailed description of this image. 
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Fig. D.2: Major axis PV-diagrams for NGC 7252E (top), NGC 7252NW (middle), and VCC 2062 (bottom). See Appendix 
D for a detailed description of this image. 
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